


NASA 

Reference 

Publication 

1275 

1992 


SAGE I Data User’s Guide 


Leonard R. McMaster 
and William P. Chu 
Langley Research Center 
Hampton, Virginia 

Michael W. Rowland 
ST Systems Corporation (STX) 
Hampton, Virginia 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Program 



Preface 

The data products from the Stratospheric Aerosol and Gas Experiment I (SAGE I), flown on 
the Applications Explorer Mission 2 (AEM-2) satellite, have been archived on magnetic tapes 
at the National Space Sciences Data Center, NASA Goddard Space Flight Center, Greenbelt, 
Maryland 20771, and are now available to researchers. The AEM-2 satellite was launched on 
February 18, 1979, but because of the failure of the spacecraft power subsystem, the mission was 
terminated on November 11, 1981. Over the operational life of the spacecraft, the instrument 
collected 33 months of scientific data, including vertical profiles of stratospheric aerosol, ozone, 
and nitrogen dioxide. This document is intended to serve as a guide for using these data sets in 
scientific investigations of stratospheric chemistry related to aerosol, ozone, nitrogen dioxide, 
dynamics, and climate change. There are brief descriptions of the instrument operation; data 
collection, processing, and validation; and some descriptions of the scientific analyses that 
have been conducted. The SAGE I data products, containing the aerosol extinction data at 
1000 and 450 nm, the ozone concentration data, and the nitrogen dioxide concentration data, 
are described in detail in appendixes A, B, and C. 

Over the years, the SAGE I experiment’s development and research have been guided by 
the SAGE I Science Team, made up of the following people: 

M. P. McCormick, Principal Investigator, NASA Langley Research Center (LaRC); R. Craig, 
Florida State University; D. M. Cunnold and G, W. Grams, Georgia Institute of Technology, 
B. M. Herman, University of Arizona; M. Hirono, Kyoto University; D. E. Miller, British 
Meteorological Office; D. G. Murcray, University of Denver; T. J. Pepin, University of Wyoming; 
W. G. Planet, National Oceanic arid Atmospheric Administration National Meteorological 
Center; and P. B. Russell, NASA Ames Research Center. 

The efforts of the LaRC SAGE 1 Data Processing Team arc also gratefully acknowledged. 
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Summary 

The Stratospheric Aerosol and Gas Experiment I 
(SAGE I) instrument was launched on board the 
Applications Explorer Mission-2 (AEM-2) satellite, 
which was dedicated to the SAGE I mission. SAGE I 
provided global data collection of aerosol vertical ex- 
tinction profiles, ozone vertical concentration pro- 
files, and nitrogen dioxide vertical concentration pro- 
files. The SAGE I instrument was a Sun photometer 
that measured the attenuation of solar radiation in 
four spectral regions through the Earth's atmosphere 
during spacecraft sunrise and sunset. The solar radi- 
ance data were combined with spacecraft ephemeris 
and National Oceanic and Atmospheric Administra- 
tion (NOAA) meteorological data and then numer- 
ically inverted to yield altitude profiles of aerosol 
extinction at wavelengths of 1000 and 450 nm and 
altitude profiles of ozone and nitrogen dioxide con- 
centration. The SAGE I aerosol data were validated 
by comparison with correlative lidar and dustsonde 
in situ measurements, the ozone data were vali- 
dated by comparison with balloon electrochemical 
cell (ECC) ozonesonde and rocket measurements, 
and the nitrogen dioxide measurements were com- 
pared with climatology. These data are currently 
archived at the National Space Sciences Data Cen- 
ter (NSSDC). This publication describes the SAGE I 
experiment, instrument characteristics, and mode of 
operation; outlines the method of the data inversion; 
and explains to potential users the format of the data, 
products archived at NSSDC. Results of the data val- 
idation and examples of data products are also pre- 
sented to demonstrate the quality of the data and 
some of its applications in atmospheric studies. 

Introduction 

SAGE I was the second in a series of satellite 
experiments using the solar occultation technique to 
monitor and study stratospheric trace constituents. 
The Stratospheric Aerosol Measurement II (SAM II), 
the first automatic instrument to use this technique 
from space, was launched October 24, 1978, on the 
Nimbus 7 spacecraft (McCormick et al. 1979). The 
SAM II instrument uses a single spectral channel 
centered at a wavelength of 1000 nm for monitoring 
stratospheric aerosols. Because of the Nimbus 7’s 
Sun-synchronous orbit, SAM II collects stratospheric 
aerosol data exclusively in the two polar regions at 
latitudes extending from 64° to 80° N and from 
64° to 80° S. At the time this paper was prepared, 
SAM II was still operational and in its 13th year of 
data collection. 

The Stratospheric Aerosol and Gas Experiment I 
(SAGE I) was launched February 18, 1979, aboard 


the Applications Explorer Mission-2 (AEM-2) satel- 
lite (McCormick et al. 1979). The SAGE I instru- 
ment had four spectral channels centered at wave- 
lengths of 1000, 600, 450, and 385 nm for measuring 
the atmospheric extinction on account of aerosols, 
ozone, and nitrogen dioxide. The AEM-2 satellite 
was placed in an orbit of approximately 600 km at 
an inclination of 56° to extend the latitudinal cov- 
erage for the solar occultation measurements from 
79° S to 79° N. The SAGE I instrument collected 
data for almost 3 years until the AEM-2 satellite 
power subsystem failed. The processed SAGE I data 
have been archived at the National Space Sciences 
Data Center (NSSDC), NASA Goddard Space Flight 
Center, Greenbelt, Maryland 20771. 

The most recent of these solar occultation ex- 
periments, the Stratospheric Aerosol and Gas Ex- 
periment II (SAGE II), uses seven spectral channels 
for measuring the atmospheric extinction of aerosols, 
ozone, nitrogen dioxide, and water vapor and is in 
an orbit similar to that of SAGE I, which allows 
data collection at latitudes from 80° S to 80° N. The 
SAGE II instrument was launched October 5, 1984, 
aboard the Earth Radiation Budget Satellite (ERBS) 
and was still collecting data after more than 7 years 
of operation at the time this paper was prepared 
(McMaster 1986). As the SAGE II data are pro- 
cessed and validated, they are also being archived at 
the NSSDC. 

AEM-2 SAGE I Mission 

The scientific objective of SAGE I was to develop 
a global stratospheric aerosol, ozone, and nitrogen 
dioxide data base that could be used for the investi- 
gation of the spatial and temporal variations of these 
species caused by seasonal and short-term meteoro- 
logical variations, atmospheric chemistry and micro- 
physics, and transient phenomena such as volcanic 
eruptions. The data base could also be used for the 
study of trends, atmospheric dynamics and trans- 
port, and potential climatic effects. The SAGE I sen- 
sor was designed to measure the attenuation of solar 
radiation resulting from atmospheric aerosol, ozone, 
and nitrogen dioxide at four spectral regions through 
the Earth's atmosphere during each spacecraft sun- 
rise and sunset (fig. 1). The SAGE I sensor was a 
dedicated instrument on the AEM-2 satellite which, 
had its orbit tailored to maximize geographic cover- 
age for solar occultation sampling. The AEM-2 satel- 
lite, with its orbital altitude of approximately 600 km 
and angle of inclination of 56°, provided 15 sunrise 
and 15 sunset measurements each day, with succes- 
sive sunrise or sunset measurements occurring about 
25° of longitude apart. The latitude of the sunrise or 


sunset measurements varied slowly from day to day 
and provided near-global coverage every 3 to 4 weeks. 
As the seasons changed, the maximum latitude of the 
SAGE I measurements ranged from 79° S to 79° N. 
Figure 2 shows the latitudinal coverage of the SAGE I 
mission over the 33 months of instrument operation. 
Because of the spacecraft power subsystem problems 
that reduced the power capacity of its batteries, sun- 
rise measurements from SAGE I were terminated 
after August 1979. 



Figure L Solar occultation measurement geometry. 

Instrument Description and Operation 

The SAGE I instrument was a four-channel Sun 
photometer that used a Cassegrainian telescope, 
holographic grating, and four silicon photodiodes to 
define the four-spectral-channel bandpass. Figure 3 
is a schematic of the SAGE I instrument. Solar ra- 
diation was reflected off a scan mirror into the tele- 
scope, with an image of the Sun being formed at the 
focal plane. The instrument’s instantaneous field of 
view, defined by the aperture on the focal plane, was 
a 30-arcsec circle that produced a vertical resolution 
at the tangent point of about 0.5 km. Radiation pass- 
ing through the aperture was transferred to the spec- 
trometer section of the instrument which contained 
the holographic grating and four separate detector 
systems. The holographic grating dispersed the in- 
coming radiation into the four spectral regions cen- 
tered at wavelengths of 1000, 600, 450, and 385 nm. 
Slits on the Rowland circle of the grating defined 
the spectral bandpass of the four spectral channels. 
The bandpasses were 50, 30, 20, and 30 nm, respec- 
tively, for the above-mentioned wavelengths. The en- 
tire imaging and spectrometer system was inside the 
azimuth gimbal to allow the instrument to be pointed 
at the Sun without image rotation. The azimuth gim- 
bal could be rotated over 360° so that measurements 
could be made at any azimuthal angle. 
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(a) 1979. 



Month 


(b) 1980. 
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(c) 1981. 

Figure 2. Latitudinal coverage for SAGE I measurements. 
Pluses indicate sunrise measurements and diamonds in- 
dicate sunset measurements. 
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Figure 3. Schematic of SAGE I instrument. (From 
McCormick et al. 1979.) 


The operation of the instrument, during each sun- 
rise and sunset measurement, was totally automatic. 
Prior to each sunrise or sunset, the instrument was 
rotated in azimuth to its predicted solar acquisition 
position. When the Sun entered the instrument’s 
field of view, the instrument adjusted its azimuth po- 
sition to lock onto the radiometric center of the Sun 
to within ±45 arcsec and then acquired the Sun by 
rotating its scan mirror to the proper elevation an- 
gle. As the Sun traversed between the horizon and 
the tangent height of 150 km, the elevation mirror 
scanned vertically across the solar disk. (See fig. 4.) 
The two solid lines in the figure represent the top 
and bottom edges of the solar disk during a sunset 
event. The gradual shrinking of the vertical shape of 
the solar image at the lower altitudes was caused by 
increasing atmospheric refraction. The dashed line 
represents the relative motion of the instrument ele- 
vation mirror as it scanned the Sun at a nominal rate 
of 15 arcmin/sec. The radiometric channel data were 
sampled at a rate of 64 samples per second per chan- 
nel, digitized to 12-bit resolution, and recorded for 
later transmission back to Earth. Additional SAGE I 
instrument information can be found in McCormick 
et al. (1979). 



Figure 4. Typical SAGE I elevation scanning of solar disk for 
a sunset measurement event. 

Data Processing 

The SAGE I measurement data were pro- 
cessed with the algorithm developed by Chu and 
McCormick (1979) and its updated revision as dis- 
cussed by Chu et al. (1989). In this section, a brief 
outline of the SAGE I algorithm is presented. 

The measured radiance from the SAGE I instru- 
ment can be related to the atmospheric optical prop- 
erties with the Lambert-Beer relation: 

H X (f) = J F (0, <t>) S {6, <j>, t) exp[-r A (0)] dQ (1) 


where 


H\(t) measured radiance at time t for 

channel wavelength A 

F(6 , <fi) instrumental field-of-vicw function 

*?(#, t) extraterrestrial solar radiance 

profile within the SAGE I spectral 
bandwidth at time t corrected for 
atmospheric refraction effects 


t\(6) optical thickness of atmosphere for 

elevation angle 6 and wavelength A 

< p azimuthal angle 

Q total solid angle 


Each elevation angle 9 corresponds uniquely to an 
atmospheric tangent height ht , the optical thickness 
t\(6) can be related to atmospheric extinction prop- 
erties through the following equation: 


r x {h t ) = 


/ 


[da (h) + 0 O3 (h) + P N q 2 (h) + /3 Ray (ft)] dp 


( 2 ) 


where 

(i n (h) profile of aerosol extinction versus 

altitude 
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/?0 3 (h) profile of ozone extinction versus 
altitude 

0NO 2 (h) profile of nitrogen dioxide extinction 

versus altitude 

/? R ay(/0 profile of Rayleigh extinction versus 

altitude 

dp differential path length through 

atmosphere 

The integral is evaluated from the spacecraft position 
to the Sun. 

The processing of the SAGE I data to generate 
the various species profiles was done in three steps. 
First, the measured radiance data were reduced with 
the spacecraft ephemeris, atmospheric refraction cal- 
culation, and the Sun-scan data and put into profiles 
of limb optical thickness r\{ht) as a function of tan- 
gent height ht in the atmosphere for each channel 
centered at wavelength A. The high- altitude solar- 
scan profiles for each channel were used as calibrated 
solar-limb profiles out of the atmosphere. 

The second step was to subtract the estimated 
Rayleigh contribution along each limb path for each 
channel and to separate the four spectral optical 
thickness profiles into aerosol optical thickness pro- 
files at 1000 and 450 nm, an ozone optical thick- 
ness profile at 600 nm, and a nitrogen dioxide op- 
tical thickness profile at 385 nm with the method 
described by Chu and McCormick (1979) and Chu 
et al. (1989). The Rayleigh contributions were com- 
puted from coincident temperature versus height pro- 
files that were provided by the National Oceanic and 
Atmospheric Administration (NOAA). 

The third step was to invert the species’ optical 
thickness profiles into extinction profiles. Through 
division of the atmosphere into N homogeneous lay- 
ers, the integral equation can be reduced to a system 
of linear equations as follows: 

N 

T ai = ^ " Pij (^) 

3 

where 

r a i measured limb optical thickness at ith layer 
for species a 

p^ path length of the Sun’s ray in j th layer 
with its tangent height at ith layer 

3 a j averaged extinction coefficient for species a 
in jth layer 


Equation (3) is then inverted with Twomey’s modifi- 
cation of Chahine’s nonlinear relaxation algorithm 
as described by Chu and McCormick (1979) and 
Chu et al. (1989). The SAGE I inversion algo- 
rithm generated vertical aerosol extinction profiles 
at 1000 and 450 nm and ozone and nitrogen dioxide 
vertical concentration profiles. 

An example of the SAGE I aerosol data inversion 
is shown in figure 5. The aerosol extinction profiles 
of 1000 nm that are shown are from SAGE I radiance 
data collected April 24, 1979, at 0545 and 0722 LCT. 
They illustrate the background aerosol condition ver- 
sus the presence of volcanic plume from the eruption 
of Soufriere volcano (McCormick et al. 1982). The 
error bars shown in the figures are estimates of the 
random errors associated with the calculations of the 
aerosol extinction that result from both the radio- 
metric measurement and the mathematical inversion. 
These errors (also archived at NSSDC along with the 
constituent profiles) include contributions from ran- 
dom measurement and inversion noise, NOAA tem- 
perature uncertainties, and altitude uncertainties. 
Below an altitude of 25 km, where aerosol extinction 
exceeds molecular extinction by 50 percent, the total 
error in the retrieved aerosol extinction coefficient is 
typically less than 10 percent. Therefore, even under 
most background or nonvolcanic conditions, the ex- 
tinction resulting from stratospheric aerosols can be 
measured to accuracies within 10 percent. 

An example of the SAGE I ozone data inver- 
sion is shown in figure 6. The dashed-line profile 
of ozone partial pressure is from SAGE I radiance 
data collected on April 2, 1979. It illustrates a typ- 
ical mid-latitude springtime ozone profile. The pro- 
file shown by the solid line in the figure is a coin- 
cident electrochemical cell (ECC) ozonesonde pro- 
file obtained at Garmisch-Partenkirchen, Germany 
(Reiter and McCormick 1982). The error bars asso- 
ciated with the SAGE I ozone profile are estimates 
of the random errors from the calculation of the ozone 
concentration profile. These errors are estimated 
similar to the errors in the aerosol extinction pro- 
files, Generally, the uncertainties of the ozone con- 
centration profiles retrieved from the SAGE I mea- 
surements are less than 10 percent from cloud- top 
height to a geometric altitude of about 45 km. 

An example of the SAGE I nitrogen dioxide data 
inversion with the 385-nm channel is shown in fig- 
ure 7. (See also WMO 1981.) Chu and McCormick 
(1986) suggested that measurement noise and diffi- 
culties associated with removing the Rayleigh, ozone, 
and aerosol contributions from this single wavelength 
retrieval led to an accuracy of approximately 20 per- 
cent, as shown by the error bars. The effect of the 
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(a) Background conditions. 0722 LCT; long. 85.4° W; 
lat. 19.8° N. 



(b) Perturbed conditions resulting from La Soufriere volcano 
eruption. 0545 LCT; long. 61.0° W; lat. 19.5° N. 

Figure 5. SAGE I aerosol extinction profile on April 24, 1979. 
Horizontal bars indicate 1 a error bar. (From McCormick 
et al. 1982.) 

inhomogeneous distribution of nitrogen dioxide re- 
sulting from its diurnal variation along the slant tan- 
gent path has been neglected in the retrieval of ni- 
trogen dioxide vertical profiles. This effect has been 



Figure 6. SAGE I measured ozone profile and coincident ECC 
ozonesonde profile obtained at Garmisch-Partenkirchen, 
Germany. Profiles obtained on April 2, 1979. ECC 
ozonesonde data starting at 1645 GMT; SAGE T mea- 
surements at 1743 GMT. (From Reiter and McCormick 
1982.) Horizontal bars indicate l<r error bar. 



Mixing ratio, ppbv 


Figure 7. SAGE I sunset nitrogen dioxide data and previous 
measurements of nitrogen dioxide as shown in WMO 
(1981). Horizontal bars indicate 1 a error bar. (From Chu 
and McCormick 1986.) 

analyzed by Kerr, Evans, and McConnell (1977) and 
the correction factors were found to be small. The 
SAGE II instrument uses two channels with rela- 
tively narrow bandwidths to make the nitrogen diox- 
ide measurement, and this retrieval results in a signif- 
icant improvement in the profile accuracy from that 
experiment (Cunnold et al. 1991). 
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SAGE I data were originally processed and 
archived at NSSDC in 1985, with corrections to 
NO A A temperature data at upper pressure levels as 
tabulated by Geller, Wu, and Gelman (1983). How- 
ever, the tabulated correction factors were later de- 
termined to have sign errors (Gelman et. al. 1986). 
The effects on the original SAGE I results were most 
pronounced in low-latitude-high-altitude regions. As 
part of the activities of the NASA World Meteo- 
rological Organization (WMO) Ozone Trends Panel 
in 1987 (Watson et al. 1988), SAGE I data were 
reprocessed with the corrected temperature adjust- 
ment factors. This user’s guide describes the revised 
SAGE I data set that has replaced the old SAGE I 
data set at NSSDC. 


Data Validation 

Before being archived, the SAGE I data were vali- 
dated through an extensive correlative measurements 
program. Two major correlative measurement exper- 
iments were conducted — the first over Natal, Brazil, 
in April 1979, and the second over Poker Flat Re- 
search Range, Alaska, in July 1979. In each exper- 
iment, correlative aerosol and ozone measurements 
were made nearly coincident in space and time with 
the SAGE I measurements. The correlative aerosol 
data consisted of SAM II-measured aerosol extinction 
profiles, airborne lidar-measured aerosol backscatter 
data, balloon-borne optical particle counter (dust- 
sonde) data, and other in situ particle counter mea- 
surements. These correlative aerosol measurements 
were collected and converted to aerosol extinctions 
for comparison with SAGE I aerosol data (Russell, 
McCormick, et al. 1981; Russell, Swissler et al. 1981; 
Russell et al. 1984). The correlative ozone data con- 
sisted of ECC ozonesondes measurements and rocket 
ozone data (McCormick et al. 1984). Figure 8 shows 
the aerosol data from the Poker Flat experiment, 
along with the SAGE I, the SAM II, the dustsonde 
(DS), filter, and the Ames Wire Impactor (AWI) 
data (Russell et al, 1984). Other comparisons of 
the SAGE I and SAM II aerosol extinction mea- 
surements have also been made (Yue, McCormick, 
and Chu 1984). Figure 9 shows the SAGE I ozone 
measurements and the ECC ozonesonde and rocket 
measurements at Poker Flat and Wallops Island, 
Virginia (McCormick et al. 1989). Comparisons 
of the SAGE I ozone data with other ECC data 
have also been made (Reiter and McCormick 1982). 
Because of the unavailability of any generally ac- 
cepted nitrogen dioxide monitor, the SAGE I nitro- 
gen dioxide measurements were compared with avail- 
able climatology as well as noncoincident balloon 


soundings and ground-based measurements (Chu and 
McCormick 1986). 

These comparisons between SAGE I aerosol, 
ozone, and nitrogen dioxide measurements and cor- 
relative data demonstrate that the SAGE I measure- 
ments agree with the correlative data to within their 
measurement and conversion uncertainties. These 
results clearly support the validity of the SAGE I 
aerosol, ozone, and nitrogen dioxide data. 

Data Products 

The basic archived products from the SAGE I 
data processing are the aerosol profile tape, the 
ozone profile tape, and the nitrogen dioxide pro- 
file tape. These data products are archived and 
can be requested on magnetic tape from the Na- 
tional Space Sciences Data Center (NSSDC), NASA 
Goddard Space Flight Center, Code 933, Green- 
belt, Maryland 20771. (Foreign requests should be 
sent to Code 930.2.) Each request should specify 
the experiment data desired (aerosol, ozone, or ni- 
trogen dioxide), the NSSDC identification number 
(79-013A-01C), and the time period of interest. All 
33 months of SAGE I data have been archived. ; 

The aerosol profile tape contains the altitude 
profiles (plus the corresponding error estimates) of 
aerosol extinction coefficients at wavelengths of 1000 
and 450 nm, Rayleigh extinction coefficients at wave- 
lengths of 1000 and 450 nm, and the extinction ratio 
of aerosol extinction coefficient to Rayleigh extinc- 
tion coefficient at a wavelength of 1000 nm. Auxiliary 
information that was used for processing the SAGE I 
data, such as NOAA's temperature-versus-altitude 
data, ephemcris-calculated latitude and longitude at 
the tangent location, beta angle of the spacecraft or- 
bit (angle between the spacecraft orbit plane and the 
Earth-Sun direction), and event duration, is also in- 
cluded on the profile tape. There is one aerosol profile 
tape per year that covers 12 months of data from De- 
cember of 1 year to November of the next year (winter 
to fall). The tapes are formatted as 6250-bpi, 9-track 
tapes and are archived at NSSDC. The aerosol profile 
tape record format is described in appendix A. - 

The ozone profile tapes contain the sunrise and 
sunset altitude profiles of ozone number density and 
mixing ratio, plus their corresponding error esti- 
mates. They contain the same auxiliary information 
as the aerosol profile tapes. There is one ozone pro- 
file tape per year that covers 12 months of data from I 

December of 1 year to November of the next year. I 

The tapes are formatted as 6250-bpi, 9-track tapes ■ 

and are archived at NSSDC. The ozone profile tape ! 

format is described in appendix B. j 
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(a) July 16. (b) July 17. (c) July 19. 

Figure 8. SAGE I aerosol extinction coefficient retrieval and correlative aerosol measurements collected during 1979 validation 
experiment over Poker Flat Research Range, Alaska. Arrows in figure indicate height of tropopause. A — 1.0 /im. 
DS = Dustsonde; AWI = Ames Wire Impactor. (From Russell et al. 1984.) 
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Ozone mixing ratio, ppmv 


Figure 9. SAGE I ozone data with coincident ECC ozonesonde data and rocket ozonesonde data. Horizontal bars indicate lcr 
error bar. (From McCormick et al. 1989.) 
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The nitrogen dioxide profile tapes contain SAGE I 
nitrogen dioxide data in the same format as the ozone 
profile tapes. Similarly, there is one nitrogen dioxide 
data tape per year that covers 12 months of data from 
December of 1 year to November of the next year. 
The nitrogen dioxide profile tape format is described 
in appendix C. 

Appendixes A, B, and C refer to the meteorologi- 
cal data notes of appendix D to explain the meteoro- 
logical fields supplied by the National Meteorological 
Center in the data records of the tapes. Appendix D 
also explains how the meteorological data were 
handled by the SAGE I processing team. 

Scientific Studies 

The SAGE I data have been used to perform a va- 
riety of atmospheric studies of atmospheric chemistry 
and dynamics, history of volcanic perturbation, and 
ozone and nitrogen dioxide long-term trends. The 
following discussion briefly describes some results of 
these investigations that used the SAGE I aerosol, 
ozone, and nitrogen dioxide data. 

SAGE I aerosol data provided the first global de- 
scription of volcanic perturbation of the stratospheric 
dust layer during the 33 months that SAGE I was in 
operation (Kent and McCormick 1984). Figure 10 
shows the temporal variation of the optical depth 
of global stratospheric aerosols from 1979 to 1981, 
based on the combined SAGE I and SAM II aerosol 
measurements. It shows the effects of stratospheric 
aerosol loading resulting from the different volcanic 


injections on both the Northern Hemisphere and 
the Southern Hemisphere. The two wavelengths 
of aerosol extinction coefficients obtained from the 
SAGE I measurements have been used to estimate 
the stratospheric aerosol size information, spatial 
distribution, and temporal evolution (Lenoble and 
Pruvost 1983; Yue and Decpak 1983; Yue and 
Deepak 1984). 

The primary goal of the SAGE I mission was 
to monitor stratospheric trace species. However, 
the SAGE I instrument did make measurements in 
the troposphere at the cloud-free regions. Stud- 
ies of the aerosols in the free troposphere, as mea- 
sured by SAGE I, have been performed to research 
seasonal variability (Kent et al. 1988). Also, the 
SAGE I tropospheric data can be used to gener- 
ate a high-altitude cloud climatology (Woodbury and 
McCormick 1986). Figure 11 illustrates a global map 
that was developed by the SAGE I data which shows 
the frequency of occurrence of the high-altitude 
clouds. 

SAGE I aerosol, ozone, and nitrogen dioxide data 
have also been used for studies of atmospheric dy- 
namics. For example, high-latitude aerosol, ozone, 
and nitrogen dioxide data have been used to study 
stratospheric sudden warming phenomena (Wang 
and McCormick 1985; Chu and McCormick 1986). 
SAGE I ozone data have been used to study the lat- 
itudinal transport in the atmosphere caused by wave 
activities (Wang, McCormick, and Chu 1983). 
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Figure 10. Temporal variation of optical depth of mean hemispheric stratospheric aerosols based on combined SAGE I and SAM II 
measurements from 1979 to 1981. Mass loading conversion factor = 1.10 x 10^ in^ kg b (From Kent and McCormick 1984.) 
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Figure 11. Average frequency of cirrus cloud occurrence during full 33 months of SAGE I operation. (From Woodbury and 
McCormick 1986.) 


SAGE I ozone and nitrogen dioxide data have also 
been used to study atmospheric chemistry anoma- 
lies related to nitrogen dioxide photochemistry (Chu 
and McCormick 1986) and ozone chemistry (Wang, 
McCormick, and Chu 1983). 

The SAGE I ozone data have been compared with 
SAGE II ozone data to determine the global ozone 
trends from 1979 to late 1980 (Watson et al. 1988; 
McCormick, Veiga, and Zawodny 1989). It was this 
comparison between SAGE I and SAGE II ozone 
data that demonstrated stratospheric ozone at 50 km 
did not decrease significantly over the period stud- 
ied contrary to the trend results based on other 
satellite observations (Herman, Hudson, and Ser- 
afino 1990; Chandra et al. 1990). The SAGE I and 
SAGE II ozone data have been regarded as the most 
accurate ozone profile data measured from space 
(WMO 1990). 


Concluding Remarks 

This report has described the Stratospheric Aero- 
sol and Gas Experiment I (SAGE I) instrument; 
its data collection, processing, and validation tech- 
niques; and the aerosol, ozone, and nitrogen diox- 
ide archival products used in scientific investigations 
of various atmospheric sciences. These data prod- 
ucts are archived and can be requested on magnetic 
tape from the National Space Sciences Data Cen- 
ter (NSSDC), NASA Goddard Space Flight Cen- 
ter, Code 933, Greenbelt, Maryland 20771. All 
33 months of SAGE I data have been archived. 


NASA Langley Research Center 
Hampton, VA 23681-0001 
May 28, 1992 
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Appendix A 

SAGE I Aerosol Profile Record Format 
Record Format 


CYBER a 

words 

(60-bit) 

Size 

Field content description 

Note 

(b) 

40-km reference data 

0001 


L 

Event Date (yymmdd.0) 

1 

0002 



Event Time (hhmmss.0) 

1 

0003 



Subtangent Latitude (0.0° ± 90.0°) 


0004 



Subtangent Longitude (0.0° ± 180.0°) 


0005 



Spacecraft- Referenced Event Type (0.0 = Sunrise; 1.0 = Sunset) 

2 

0006 



Earth- Referenced Event Type (0.0 = Sunrise; 1.0 = Sunset) 

2 

0007 



Spacecraft Beta Angle (0.0° ±61.0°) 

3 

0008 



Coded Time of Year (ddd.fract) 

4 

NMC meteorological data (see appendix D) 

0009 0033 

25 

Temperature, K 


0034-0058 



Temperature Error, K 


00590083 



Geometric Altitude, m 


0084 0108 



Air Density, g/m 3 


0109 0133 

' 


Air Density Error, percent 


0134 

1 

Temperature Correction Value for 5.0-mbar Level, K 


0135 



Temperature Correction Value for 2.0-mbar Level, K 


0136 



Temperature Correction Value for 1.0-mbar Level, K 


0137 



Temperature Correction Value for 0.4-mbar Level, K 


0138 



“Meteorological Data Not Complete” Flag (0 = Complete; 1 — Incomplete) 


0139 



“Start of Model Meteorological Data” Array Index Pointer (1 19) 


0140 



Model Meteorological Data Selection Code (ssll) 


0141 



Revision Date of LaRC Meteorological Model (yymmdd.0) 


NASA LaRC processing information 

0142 

1 

LaRC Driver Revision Level 


0143 



LaRC Transmission Revision Level 


0144 



LaRC Inversion Revision Level 


0145 



LaRC Event Tag (yymmddhhmm.sq) 

5 

0146 



LaRC Processing Date (yymmdd.0) 

1 

0147 



LaRC Processing Time (hhmmss.0) 

1 

0148 



Mean Subtangent Altitude for Event Limb Calibration, km 

6 

0149 

' 


Value Designated as the Data Fill Number for this Event 

7 

Event ground-track slew data 

0150 0157 

8 

Subtangent Altitude, km 


0158 0165 

8 

Corresponding Latitude (0.0° ± 90.0°) 


0166 0173 

8 

Corresponding Longitude (0.0° ± 180.0°) 


0174 

1 

Time Span of Data from Levels 1 to 70, sec 


Altitude and meteorological data for profile arrays 

0175 0244 

70 

Geometric Altitude, km 


0245-0314 

70 

Corresponding Pressure, mbar 


0315 0384 

70 

Corresponding Temperature, K 


0385 0390 

6 

Spare 



“Trademark of Control Data Corporation. 
6 Data field notes given after table. 
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CYBER 0 

words 

(60-bit) 

Size 

Field content description 

Note 

(b) 

Quality estimations of channel optical depth profile 

0391 

1 

1000-nm- Wavelength Quality Factor 

8 

0392 



Spare 


0393 



Spare 


0394 



Spare 


0395 



450-nm- Wavelength Quality Factor 

8 

0396 



Spare 


0397 



Spare 


0398-0400 

3 

Spare 


Rayleigh extinction profiles 

0401 0460 

60 

1000-nm Rayleigh Extinction, km - 


0461-0520 



1000-nm Rayleigh Extinction Error, km -1 


0521-0580 



Spare 


0581-0640 



Spare 


0641 0700 



450-nm Rayleigh Extinction, km -1 


0701 0760 



450-nm Rayleigh Extinction Error, km -1 


0761-0820 



Spare 


0821-0880 

- 


Spare 


Aerosol profiles 

0881 0940 

60 

1000-nm Extinction, km -i 


0941 1000 



1000-nm Extinction Error, km -1 


1001 1060 



Spare 


1061-1120 



Spare 


1121-1180 



450-nm Extinction, km” 1 


1181-1240 



450-nm Extinction, km -1 


1241-1300 



Spare 


1301-1360 



Spare 


1361-1420 



1000-nm Extinction Ratio 

9 

1421-1480 


«• 

1000-nm Extinction Ratio Error 

9 

1481-1488 

8 

Spare 


End of event record I 


° Trademark of Control Data Corporation. 
& Data field notes given after table. 
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Record Format Notes 

General notes: 

• Each field of the event record contains one 60-bit CYBER a floating point number. 

• All time and data references are to GMT, except fields 146 and 147, which are LaRC processing time. 

• All latitudes and longitudes are given at the event subtangent point. 

• If any field in the event record is considered invalid, or has missing data, a fill value will be placed 
in that field. For each event record, that fill value can be found in field 149. (See data field note 7 
below.) 

• Each profile level bin is centered at the 0.5-km point within a range of 1.0 km. 

Data field notes: 

1. The “yymmdd.O” and “hhmmss.O” fields are generated by the respective FORTRAN statements: 

DATE = FLOAT (IYY* 10000 + IMM*100 + IDD) 

TIME = FLOAT (IHH* 10000 + IMM*100 + ISS) 

2. Spacecraft- Referenced Event Type and Earth- Referenced Event Type fields are normally the same 
type, but if the absolute value of the spacecraft beta angle is close to 61°, their values may be opposite. 
The Earth-Referenced Event Type field is based on a ground-observer’s viewpoint. 

3. The Spacecraft Beta Angle field is defined as the angle generated by the intersection of the Earth-Sun 
vector and the spacecraft orbit plane. 

4. The Coded Time of Year field is the time at the beginning of the event (not the same time as for 
fields 1 and 2) and is generated by the FORTRAN statement: 

CODTIME = FLOAT(DOY) + (SOD/86400.0) 
where 

DOY = Day of year (1-366) 

SOD = Seconds of the day (0.0 86399.99...) 

5. The LaRC Event Tag field is generated with statements similar to those in note 1. The “.sq” at the 
end of the value is the event number of the day divided by 100. 

6. The Mean Subtangent Altitude for Event Limb Calibration field contains the altitude at which data 
for the exoatmospheric solar image were gathered for use in solar limb normalization for the event. 

7. The Value Designated as the Data Fill Number for This Event field must be used to determine what 
data in the event record are valid. If any field other than this one contains this number, that field has 
no valid information and should not be used by the investigator. 

8. The Quality Factor fields for each wavelength are equal to 1.0 minus the summation of the optical 
depth errors at each profile level from 20.5 to 59.5 km. In cases where a 40-km span cannot be realized, 
the quality factor is proportioned to a 40-km span to allow a better comparison across wavelengths 
and other events. 

9. Extinction ratio = (Aerosol extinction + Rayleigh extinction) /Rayleigh extinction. 
a Trademark of Control Data Corporation. 
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Appendix B 

SAGE I Ozone Profile Record Format 

Record Format 


CYBER 0 

words 

(60-bit) 

Size 

Field content description 

Note 

(b) 

40-km reference data 

0001 

] 


Event Date (yymmdd.0) 

1 

0002 



Event Time (hhmmss.0) 

1 

0003 



Subtangent Latitude (0.0° =t 90.0°) 


0004 



Subtangent Longitude (0.0° ± 180.0°) 


0005 



Spacecraft- Referenced Event Type (0.0 = Sunrise; 1.0 = Sunset) 

2 

0006 



Earth-Referenced Event Type (0.0 = Sunrise; 1.0 = Sunset) 

2 

0007 



Spacecraft Beta Angle (0.0° ±61.0°) 

3 

0008 

* 


Coded Time of Year (ddd.fract) 

4 

NMC meteorological data (see appendix D) 

0009 0033 

25 

Temperature, K 


0034-0058 



Temperature Error, K 


0059 0083 



Geometric Altitude, m 


0084 0108 



Air Density, g/m 3 


0109-0133 

- 


Air Density Error, percent 


0134 

1 

Temperature Correction Value for 5.0-mbar Level, K 


0135 



Temperature Correction Value for 2.0-mbar Level, K 


0136 



Temperature Correction Value for 1.0-mbar Level, K 


0137 



Temperature Correction Value for 0.4-mbar Level, K 


0138 



“Meteorological Data Not Complete” Flag (0 — Complete; 1 = Incomplete) 


0139 



“Start of Model Meteorological Data” Array Index Pointer (1-19) 


0140 



Model Meteorological Data Selection Code (ssll) 


0141 



Revision Date of LaRC Meteorological Model (yymmdd.0) 


NASA LaRC processing information 

0142 

1 

LaRC Driver Revision Level 


0143 



LaRC Transmission Revision Level 


0144 



LaRC Inversion Revision Level 


0145 



LaRC Event Tag (yymmddhhmm.sq) 

5 

0146 



LaRC Processing Date (yymmdd.0) 

1 

0147 



LaRC Processing Time (hhmmss.0) 

1 

0148 



Mean Subtangent Altitude for Event Limb Calibration, km 

6 

0149 


- 

Value Designated as the Data Fill Number for this Event 

7 

Event ground- track slew data 

0150-0157 

8 

Subtangent Altitude, km 


0158 0165 

8 

Corresponding Latitude (0,0° ± 90.0°) 


0166 0173 

8 

Corresponding Longitude (0.0° ± 180.0°) 


0174 

1 

Time Span of Data from Levels 1 to 70, sec 


Altitude and meteorological data for profile arrays 

0175 0244 

70 

Geometric Altitude, km 


0245 0314 

70 

Corresponding Pressure, mbar 


0315 0384 

70 

Corresponding Temperature, K 


0385-0390 

6 

Spare 



“Trademark of Control Data Corporation, 
6 Data field notes given after table. 
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CYBER" 

words 

(60-bit) 

Size 

Field content description 

Note 

(b) 

Quality estimations of channel optical depth profile 

0391 

i 

Spare 


0392 



Spare 


0393 



600-nm-Wavelength Quality Factor 

8 

0394 



Spare 


0395 



Spare 


0396 



Spare 


0397 



Spare 


0398 0400 

3 

Spare 


[ Ozone profiles j 

0401-0470 

70 

Number Density, molecules/cm 3 


0471 0540 



Number Density Error, molecules/cm 3 


0541-0610 



Volumetric Mixing Ratio, v/v 


0611-0680 

J 


Volumetric Mixing Ratio Error, v/v 


0681 0688 

8 

Spare 


End of event record 


a Trademark of Control Data Corporation. 
fe Data field notes given after table. 



Record Format Notes 


General notes: 

• Each field of the event record contains one 60-bit CYBER 0 floating point number. 

• All time and data references are to GMT, except fields 146 and 147, which are LaRC processing time. 

• All latitudes and longitudes are given at the event subtangent point. 

• If any field in the event record is considered invalid, or has missing data, a fill value will be placed 
in that field. For each event record, that fill value can be found in field 149. (See data field note 7 
below.) 

• Each profile level bin is centered at the 0.5-km point within a range of 1.0 km. 

Data field notes: 

1. The “yymmdd.O” and “hhmmss.O” fields are generated by the respective FORTRAN statements: 

DATE = FLOAT (IYY* 10000 + IMNT100 + IDD) 

TIME = FLOAT (IHH* 10000 + IMM’100 + ISS) 

2. Spacecraft-Referenced Event Type and Earth-Referenced Event Type fields are normally the same 
type, but if the absolute value of the spacecraft beta angle is close to 61°, their values may be opposite. 
The Earth- Referenced Event Type field is based on a ground-observer’s viewpoint. 

3. The Spacecraft Beta Angle field is defined as the angle generated by the intersection of the Earth-Sun 
vector and the spacecraft orbit plane. 

4. The Coded Time of Year field is the time at the beginning of the event (not the same time as for 
fields 1 and 2) and is generated by the FORTRAN statement: 

CODTIME = FLOAT(DOY) + (SOD/86400.0) 
where 

DOY = Day of year (1-366) 

SOD = Seconds of the day (0.0-86399.99...) 

5. The LaRC Event Tag field is generated with statements similar to those in note 1. The “.sq” at the 
end of the value is the event number of the day divided by 100. 

6. The Mean Subtangent Altitude for Event Limb Calibration field contains the altitude at which data 
for the exoatmospheric solar image were gathered for use in solar limb normalization for the event. 

7. The Value Designated as the Data Fill Number for This Event field must be used to determine what 
data in the event record are valid. If any field other than this one contains this number, that field has 
no valid information and should not be used by the investigator. 

8. The Quality Factor fields for each wavelength are equal to 1.0 minus the summation of the optical 

depth errors at each profile level from 20.5 to 59.5 km. In cases where a 40-km span cannot be realized, 
the quality factor is proportioned to a 40-km span to allow a better comparison across wavelengths 
and other events. 

"Trademark of Control Data Corporation. 
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Appendix C 

SAGE I Nitrogen Dioxide Profile Record Format 
Record Format 


CYBER a 

words 

(60-bit) 

Size 

Field content description 

Note 

(b) 

40-km reference data 

0001 



Event Date (yymmdd.0) 

1 

0002 



Event Time (hhmmss.0) 

1 

0003 



Subtangent Latitude (0.0° ± 90.0°) 


0004 



Subtangent Longitude (0.0° ± 180.0°) 


0005 



Spacecraft-Referenced Event Type (0.0 = Sunrise; 1.0 = Sunset) 

2 

0006 



Earth-Referenced Event Type (0.0 = Sunrise; 1.0 = Sunset) 

2 

0007 



Spacecraft Beta Angle (0.0° =t 61.0°) 

3 

0008 

' 


Coded Time of Year (ddd.fract) 

4 

NMC meteorological data (see appendix D) 

0009 0033 

25 

Temperature, K 


0034 0058 



Temperature Error, K 


0059 0083 



Geometric Altitude, m 


0084^)108 



Air Density, g/m 3 


0109 0133 

> 


Air Density Error, percent 


0134 

1 

Temperature Correction Value for 5.0-mbar Level, K 


0135 



Temperature Correction Value for 2.0-mbar Level, K 


0136 



Temperature Correction Value for 1.0-mbar Level, K 


0137 



Temperature Correction Value for 0.4-mbar Level, K 


0138 



“Meteorological Data Not Complete” Flag (0 = Complete; 1 = Incomplete) 


0139 



“Start of Model Meteorological Data” Array Index Pointer (1-19) 


0140 



Model Meteorological Data Selection Code (ssll) 


0141 



Revision Date of LaRC Meteorological Model (yymmdd.0) 


NASA LaRC processing information 

0142 


LaRC Driver Revision Level 


0143 



LaRC Transmission Revision Level 


0144 



LaRC Inversion Revision Level 


0145 



LaRC Event Tag (yymmddhhmm.sq) 

5 

0146 



LaRC Processing Date (yymmdd.0) 

1 

0147 



LaRC Processing Time (hhmmss.0) 

1 

0148 



Mean Subtangent Altitude for Event Limb Calibration, km 

6 

0149 



Value Designated as the Data Fill Number for this Event 

7 

Event ground- track slew data 

0150 0157 

8 

Subtangent Altitude, km 


0158 0165 

8 

Corresponding latitude (0.0° ± 90.0°) 


0166-0173 

8 

Corresponding Longitude (0.0° ± 180.0°) 


0174 

1 

Time Span of Data from Levels 1 to 70, sec 


Altitude and meteorological data for profile arrays 

0175 0244 

70 

Geometric Altitude, km 


0245-0314 

70 

Corresponding Pressure, mbar 


0315-0384 

70 

Corresponding Temperature, K 


0385-0390 

6 

Spare 



fl Trademark of Control Data Corporation. 
& Data field notes given after table. 
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CYBER a 

words 

(60-bit) 

Size 

Field content description 

Note 

(b) _| 

Quality estimations of channel optical depth profile 


0391 

1 

Spare 


0392 



Spare 


0393 



Spare 


0394 



Spare 


0395 



450-nm- Wavelength Quality Factor 

8 

0396 



Spare 


0397 



385-nm- Wavelength Quality Factor 

8 

0398-0400 

3 

Spare 


Nitrogen dioxide profiles 

0401-0460 

6 

0 

Number Density, molecules /crcr 


0461 0520 



Number Density Error, molecules/cm 3 


0521-0580 



Volumetric Mixing Ratio, v/v 


0581-0640 


' 

Volumetric Mixing Ratio Error, v/v 


End of event record 


“Trademark of Control Data Corporation. 
6 Data field notes given after table. 
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Record Format Notes 

General notes: 

• Each field of the event record contains one 60-bit CYBER a floating point number. 

• All time and data references are to GMT, except fields 146 and 147, which are LaRC processing time. 

• All latitudes and longitudes are given at the event subtangent point. 

• If any field in the event record is considered invalid, or has missing data, a fill value will be placed 
in that field. For each event record, that fill value can be found in field 149. (See data field note 7 
below.) 

• Each profile level bin is centered at the 0.5-km point within a range of 1.0 km. 

Data field notes: 

1. The “yymmdd.O” and “hhmmss.O” fields arc generated by the respective FORTRAN statements: 

DATE = FLOAT (IYY* 10000 -F IMM*100 + IDD) 

TIME = FLOAT (IHH* 10000 4- IMM*100 + ISS) 

2. Spacecraft-Referenced Event Type and Earth- Referenced Event Type fields are normally the same 
type, but if the absolute value of the spacecraft beta angle is close to 61° , their values may be opposite. 
The Earth-Referenced Event Type field is based on a ground-observer’s viewpoint. 

3. The Spacecraft Beta Angle field is defined as the angle generated by the intersection of the Earth-Sun 
vector and the spacecraft orbit plane. 

4. The Coded Time of Year field is the time at the beginning of the event (not the same time as for 
fields 1 and 2) and is generated by the FORTRAN statement: 

CODTIME = FLOAT(DOY) + (SOD/86400.0) 
where 

DOY = Day of year (1-366) 

SOD = Seconds of the day (0.0 86399.99...) 

5. The LaRC Event Tag field is generated with statements similar to those in note 1. The “.sq” at the 
end of the value is the event number of the day divided by 100. 

6. The Mean Subtangent Altitude for Event Limb Calibration field contains the altitude at which data 
for the exoat mospheric solar image were gathered for use in solar limb normalization for the event. 

7. The Value Designated as the Data Fill Number for This Event field must be used to determine what 
data in the event record are valid. If any field other than this one contains this number, that field has 
no valid information and should not be used by the investigator. 

8. The Quality Factor fields for each wavelength are equal to 1.0 minus the summation of the optical 
depth errors at each profile level from 20.5 to 59.5 km. In cases where a 40-km span cannot be realized, 
the quality factor is proportioned to a 40-km span to allow a better comparison across wavelengths 
and other events. 

a Trademark of Control Data Corporation. 
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Appendix D 

Meteorological Data Notes 

Meteorological data are supplied to NASA Langley Research Center (LaRC) by NO A A National Weather 
Service, Climate Analysis Branch, Washington, D.C. Data for temperature, temperature error, geometric 
altitude, air density, and air density error are provided for 18 pressure levels (1000 to 0.4 mbar) and for the 
derived tropopause pressure. The pressure levels (in mbar) correspond to the 25-element meteorological data 
arrays (1 to 25) as follows: 1000, 850, 700, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 10, 5, 2, 1, 0.4, 0.04, 0.01, 
spare, spare, spare, spare, and derived tropopause pressure. Elements 19 and 20 (0.04 and 0.01 mbar) contain 
LaRC model data for temperature and altitude only. Elements 21 to 24 (spares) contain fill values, and 
element 25 contains the NOAA-supplied tropopause information. 

If NOAA cannot supply meteorological data as above, LaRC determines the highest pressure level for which 
data are supplied and then inserts model data from the next level up to the lowest pressure level of 0.01 mbar. 
Only temperature and altitude information are supplied from using these model data. Temperature error, 
density, and density error will be filled for the corresponding levels that contain the LaRC-supplied data. 

Meteorological correction factors for temperatures at 5, 2, 1, and 0.4 mbar are already added to the value 
of the temperatures in elements 15 through 18 of the temperature array. These correction factors are only 
included in the NOAA-supplied data. If model data are in these locations, no correction factors are applied. 
The correction factors are listed in fields 134 to 137 of. the record. 

Meteorological data (fields 245 to 384 of the event record) are interpolated from the meteorological data in 
fields 9 to 133 of the record. Altitude data in fields 175 to 244 increment by 1 km with the center of each level 
at the 0.5-km point of the level bin. 

Other meteorological data information is contained in the following record locations: 

Field 0138: 0 if NOAA-supplied data are complete, 1 if incomplete. 

Field 0139: the model pointer is the array index that points to the start of LaRC-supplied model data in 
the temperature and altitude arrays of the meteorological data. 

Field 0140: model selection code (ssll) where ss is from 01 to 04 for spring to winter, and 11 is from 0 to 80 
in 10-deg increments for absolute latitude. 

Field 0141: date of revision of the LaRC-supplied model. 
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